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We report the detection of a fully-resolved, Compton-scattered emission line 

in the X-ray spectrum of the massive binary GX 301—2 obtained with the High 

Energy Transmission Grating Spectrometer onboard the Chandra X-ray Obser- 

Tj- ! vatory. The iron Ka fluorescence line complex observed in this system consists 

of an intense narrow component centered at an energy of E = 6.40 keV and a 

redward shoulder that extends down to ~ 6.24 keV, which corresponds to an 

energy shift of a Compton back-scattered iron Ka photon. From detailed Monte 

Carlo simulations and comparisons with the observed spectra, we are able to di- 

rectly constrain the physical properties of the scattering medium, including the 

electron temperature and column density, as well as an estimate for the metal 

"Sa '. abundance. 
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Introduction 



When a high energy photon propagating through a low-temperature (< 10 5 K) medium 
undergoes Compton scattering with the constituent electrons, the energy of the photon is 
modified in a way that depends on the scattering angle and the electron velocity distribution. 
In the limit where the electrons are at rest, a fraction of the photon energy is transferred to 
the electron according to the Compton formula, 



where E is the energy of the incoming photon, E\ is the energy of the outgoing photon, 9 
is the angle between the incoming and outgoing photons, and m e c 2 is the electron rest-mass 
energy (= 511 keV). The maximum energy shift per scattering due to the electron recoil is, 
therefore, AE max = 2E 2 /(m e c 2 + 2E ) for photons that are back-scattered (9 = 180°). An 
X-ray emission line propagating through a medium with substantial Compton optical depth 
( r Compton > 0.1) has a non-negligible probability of interacting with an electron, resulting in 
down-scattering of photons and, hence, producing a discernible "Compton shoulder" between 
E and E - A£ max . 

X-ray emission lines at higher energies are ideal for studying the properties of the 
Compton shoulders, since the Compton scattering opacity relative to that of photoelec- 
tric absorption is larger for higher energy photons. The iron Ka fluorescent line complex at 
E = 6.40 keV, therefore, is particularly promising, and can be produced over an extremely 
wide range in column density, which makes it ideal for diagnosing the physical properties of 
a cold medium irradiated by X-rays (Hatchett & Weaver 1977; Basko 1978; Illarionov et al. 
1979; George & Fabian 1991; Sunyaev & Churazov 1996; Matt 2002). The energy shift of an 
iron line photon due to a single Compton scattering Ai? max = 156 eV (see eq.(l) ) is, however, 
at best, comparable to the spectral resolving power of past X-ray detectors in orbit. Previ- 
ous detections of the Compton shoulder are unresolved (see, e.g., Molendi, Bianchi & Matt 
2003) and, in many cases, they are statistically marginal as well (Iwasawa, Fabian & Matt 
1997; Matt & Guainazzi 2003). The high spectral resolution of 33 eV (FWHM) at 6.4 keV 
attained by the High Energy Transmission Grating Spectrometer (HETGS) (Canizares et 
al. 2000) onboard the Chandra X-ray Observatory (Weisskopf et al. 2000) enables us, for the 
first time, to resolve the Compton shoulder of the iron Ka line (see, e.g., Kaspi et al. 2002; 
Bianchi et al. 2002). 
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The X-ray binary pulsar GX 301—2 is one of the most promising targets for studying 
the Compton shoulder associated with the 6.4 keV iron Ka line. The system consists of an 
accreting magnetized neutron star in a highly eccentric orbit (e = 0.46; Sato et al. 1986; 
Koh et al. 1997), embedded in the stellar wind from a B2 super-giant companion star (Bord 
et al. 1976). X-ray continuum photons from the neutron star ionize the K-shell electrons 
of iron in the wind, a fraction of which is followed by fluorescence emission of the 6.4 keV 
line. Previous observations of the system have shown an intense iron line with an equivalent 
width as large as ~ 1 keV (White & Swank 1984; Leahy et al. 1989; Tashiro et al. 1991; 
Endo et al. 2002). Given the large mass- loss rate of the companion (Parkes et al. 1980), a 
substantial fraction of the 6.4 keV photons can subsequently undergo Compton scattering 
prior to escaping the wind. 

In this Letter, we present spectral analysis results of the properties of the iron Ka 
line and its Compton shoulder observed in GX 301—2 with the Chandra HETGS. As we 
show, the detection of the fully-resolved Compton shoulder provides a unique opportunity 
to investigate in detail the physical state and geometrical distribution of cold (< 10 5 K) 
material surrounding an X-ray source. In § 2, we describe the observation and the procedures 
adopted for data reduction. Some basic spectral measurements and observed quantities are 
also discussed. In § 3, we describe our Monte Carlo simulator used for the spectral analyses. 
Finally in § 4, we present our measurements and discuss the implications of the results. 

2. Observation and Data Reduction 

Chandra observed GX 301—2 at three different orbital phases: (1) = 0.167-0.179, (2) 
= 0.480-0.497, and (3) = 0.970-0.982, hereafter referred to as intermediate (IM), near- 
apastron (NA), and pre-periastron (PP) phases (Pravdo et al. 1995), respectively. We focus 
primarily on the spectrum of the PP-phase, during which the X-ray luminosity, absorbing 
column density, and the iron line equivalent width are the highest (Endo et al. 2002). The 
observation of the PP-phase was performed on 2002 February 3 12:34:10 UT and continued 
for ~ 40 ksec. All of the data were processed using CIAO v2.2.1. Since the zeroth order 
image was severely piled-up, the centroid was determined by finding the intersection of the 
streak events and the dispersed events. We apply a standard spatial filter and an order 
sorting mask to extract only the first order events. The exposure map and effective area 
were calculated based on the extraction region described above. Since we are interested 
mainly in the properties of the iron line, we use data only from the High-Energy Grating 
(HEG) for the analysis, since it has a higher resolving power (AE F whm = 33 eV at 6.4 keV) 
and larger effective area (20 cm 2 ) in the iron K line region. 
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The X-ray light curve in the 1.0-10.0 keV band extracted from the HEG events is shown 
in the top panel of Fig. 1. An X-ray outburst is seen in the first half of the observation. 
Hence, we have chosen to divide the data into the first and second halves as indicated by 
the vertical dotted line, and have extracted spectra from each data segment separately. The 
observed 2-10 keV fluxes are 13.2 x 10~ 10 erg cm~ 2 s _1 and 9.0 x 10 -10 erg cm -2 s _1 
during the first and second halves, respectively. A blow-up of the iron Ka line spectra are 
shown in the bottom panels of Fig. 1. In addition to the intense iron Ka line at 6.4 keV, 
a shoulder extending towards the low energy side down to ~ 6.24 keV is clearly seen in the 
data. The equivalent widths of the iron Ka lines including the shoulders are 643 ± 20 eV 
and 486 ± 18 eV for the first and second halves, respectively. The width of the shoulder 
(AE ~ 160 eV) precisely matches what is predicted by eq. (1), which strongly indicates that 
the feature is formed primarily through single Compton scattering of the iron Ka photons. 
Changes in the shape of the Compton shoulder as well as in the shoulder flux relative to the 
un-scattered line flux, can be seen between the spectra of the first and second halves. 



3. Monte Carlo Simulation 

The Compton shoulder can be used to infer various physical parameters that characterize 
the scattering medium. The flux ratio of the shoulder to the line is determined by the metal 
abundance and the optical thickness of the scattering cloud. Its energy distribution, on the 
other hand, is sensitive to the temperature and the geometrical distribution of the scattering 
electrons. A discernible change in the profile shown in Fig. 1 implies that these physical 
parameters are variable between the first and second halves of the observation. 

In order to obtain some quantitative information from the spectra, we have constructed a 
Monte Carlo simulator to compute the emergent spectrum from an X-ray source surrounded 
by a cloud. We assume a spherical distribution of material motivated by the observed cor- 
relation between the iron emission line strength and its corresponding edge depth during 
the IM, NA, and PP phases. This assumption also significantly simplifies the simulation. 
The cloud consists of H, He and astrophysically abundant metals (C, N, O, Ne, Na, Mg, Al, 
Si, S, CI, Ar, Ca, Cr, Fe and Ni), and all of the metal abundances (Z > 2) are allowed to 
vary together relative to the cosmic values of Feldman (1992). We account for photoelectric 
absorption and subsequent fluorescent emission, as well as Compton scattering by free elec- 
tron. The angular-dependence of the Compton scattering cross section is fully accounted 
for and the electrons are assumed to have a Maxwellian energy distribution. The photons 
may suffer multiple interactions and are traced until they completely escape the cloud. The 
energy distribution of the emergent photons are then histogrammed to produce a spectrum. 
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Fig. 2 shows some of the results from the simulations for the iron line and its Compton 
shoulder with varying hydrogen column density (N H ) and electron temperature (kT c ). The 
original iron Ka± and Ka 2 photons are assumed to be distributed according to the K-shell 
photoionization rate at each radius from the central continuum source. The distribution 
of the original iron K photons is also calculated with the simulator, in which a power-law 
X-ray source with a photon index of 1.0 is assumed. In the upper panels of Fig. 2, one 
can see an increase in the scattered flux relative to the narrow line flux as iV H is increased. 
The lower panels of Fig. 2 show the temperature dependence of the shape of the Compton 
shoulder. More smearing is seen at higher kT e . Photons between 6.24 keV and 6.40 keV are 
due primarily to single-scattered photons, and the component below 6.24 keV result from 
multiple-scattering, which are important even at moderate optical depths. 

4. Discussion 

We performed spectral fits based on the simulation for the iron line and the Compton 
shoulder. A fits file of "table models" were generated from the results of the simulations, 
which were then incorporated into XSPEC vll.2. The parameters of the model are A%, 
kT e , and the metal abundances. The radial dependence of the intrinsic Ka line emissivity 
is fixed to what one expects from a photon index of T = 1.0. We have confirmed that this 
assumption produces at most a 5 % error in the emergent line profile for T between 0.0 and 
2.0. 



4.1. Metal Abundance 

We first attempt to determine the metal abundances from the total PP-phase spectrum. 
In doing this, we note that two independent constraints between the abundance and the 
hydrogen column density are available, as described below. The first one arises from the 
observed intensity ratio of the shoulder to the narrow line. For a given ratio, an increase 
in column density must be accompanied by an increase in abundance and, hence, the two 
parameters are correlated as shown in Fig. 3. The other constraint originates from the 
observed equivalent width of the line. In this case, the abundance is anti-correlated with the 
column density. From these two constraints (see Fig. 3), the metal abundance is determined 
to be 0.65-0.82 times of the cosmic value (Feldman 1992), assuming T = 1.0. The slope, 
however, is not well-measured due to the limited bandpass of the Chandra gratings. If we 
allow a range in T of 1.0-1.5, which has been seen in previous hard X-ray observations 
(Pravdo et al. 1995; Orlandini et al. 2000), the metal abundance is determined to be 0.65- 
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0.90 times of the cosmic. Note that this is consistent with that of typical OB-stars (Daflon 
et al. 2001). 

4.2. Hydrogen Column Density and Electron Temperature 

Adopting the metal abundance of 0.75 cosmic derived in § 4.1, we find the values for N-R. 
The derived values are listed in Table 1, and the models are shown by the lines superimposed 
on the data in Fig. 1. The difference in the observed Compton profile can be described by 
a change in the column density, which results in a variation in the number-of-scatterings 
distribution even at these moderate optical depths. We have also obtained upper-limits 
(90 % confidence levels) to kT c of < 3.4 eV and < 0.6 eV for the first- and second-half 
spectra, respectively. Interestingly, the derived electron temperatures are consistent with 
that expected for the stellar wind being cooled via adiabatic expansion from the surface of the 
companion star with a photospheric temperature of ~ 2 eV (Bord et al. 1976). Interestingly, 
the column density as inferred from the Compton profile, in fact, reproduces the spectrum 
in the entire HEG bandpass for a power-law photon index of 1.0. Fig. 4 shows the observed 
spectra overlaid with the simulated spectra. Reduced chi-squared values of 1.30 and 0.95 
were obtained for the first and second halves spectra, respectively. Though some residual 
flux still remain, the simulations, which are based on parameters derived from the line and 
Compton shoulder, provide fairly good descriptions of the broadband data. Assuming a 
distance of 1.8 kpc (Parkes et al. 1980), the absorption-corrected X-ray luminosities in the 
2-10 keV are 3.5 x 10 36 erg s" 1 and 1.6 x 10 36 erg s^ 1 for the first and second halves, 
respectively. 

4.3. Future Prospects 

We have demonstrated through our quantitative analysis of the Chandra HETGS data 
of GX 301—2 that the Compton shoulder provides a sensitive diagnostic for the physical 
conditions of cold material irradiated by X-ray photons. The substantial improvement in 
the spectral resolution (6 eV at 6 keV) and the effective area (~ 300 cm 2 at 6 keV) accessible 
with the X-ray Spectrometer (XRS) onboard the Astro-E2 observatory, which is planned to 
be deployed in 2005, will allow us to study the Compton shoulder profile in various classes of 
X-ray sources, such as a molecular cloud irradiated by a proto-star, an accretion disk swirling 
around a compact object, and a dusty torus surrounding an active galactic nucleus. More 
detailed geometrical information can be obtained by studying the energy distribution and, 
hence, the scattering angle distribution in the circumsource medium. Compton scattering 
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by electrons in bound systems produces spectral signatures that are distinct from those 
produced by free electrons (Sunyaev & Churazov 1996; Vainshtein et al. 1998; Sunyaev, 
Uskov & Churazov 1999), which may be detectable with future X-ray spectrometers. 

SW is supported by research fellowships of the Japan Society for the Promotion of 
Science for young Scientists. MS was supported by NASA through the Chandra Fellowship 
Program. 
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Table 1. Derived Parameters from Spectral fits of the Compton shoulder 





N u (10 23 cm- 2 ) 


T"Compton 


kT e (eV) (upper limit) 


X 2 / d.o.f. a 


First 
Second 


12.0 +3 J 


0.96^ 

o.681S:1S 


0.5 (< 3.4) 
0.0 (< 0.6) 


65.6 / 71 
82.4 / 71 



Note. — Errors and upper limits designate 90 % confidence levels. 
a Only the 6.0 - 6.6 keV region has been used in the spectral fit. 
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| : (a) Light Curve : 




6,0 6.3 6.4 6.6 6.0 6.2 6.4 6.6 
Energy[keV] Energy[keV] 

Fig. 1. — (a): The light curve of GX 301—2 in the pre-periastron (PP) phase observed with 
the Chandra HETGS in the band 1-10 keV (± 1 order of HEG). (b),(c): The spectra of the 
iron Ka region (6.0-6.6 keV) for the first and the second halves of the observation as defined 
by the dotted vertical line shown in panel (a). Each spectrum shows a clear drop in flux 
near 6.24 keV. Superimposed on the data in lines are the best-fit Monte Carlo models (see 
text and Fig. 2 ). 
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6.0 6.1 6.3 6.3 6.4 6.0 6.1 6.3 6.3 6.4 6.0 6.1 6.3 6.3 6 4 6.0 6.1 6.3 6,3 6.4 
Energy T_keV] Energy [keV] Energy [keV] Energy [keVj 



Fig. 2. — Dependence of the iron Ka line profile on the hydrogen column density (Nu) and 
the electron temperature (kT e ). The upper panels show the variation of the iron Ka line and 
its shoulder as a function iVfi for a fixed kT c at eV. The lower panels show the variation 
as a function of kT e between eV and 5 eV for a fixed iV H at 1 x 10 24 cm -2 . In these 
simulations, the metal abundances were assumed to be 0.75 times cosmic (Feldman 1992). 
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Fig. 3. — Confidence contour for 7V H vs. the metal abundance during the first half (left) 
and the second half (right). Each contour represents 68 %, 90 % and 99 % confidence level. 
The region between each pair of dashed lines are the values allowed by the observed line 
equivalent widths of 643 ± 20 eV (left) and 486 ± 18 eV (right) for two assumed values for 
the photon index (r = 1.0 and T = 1.5). The metal abundance is determined to be 0.65-0.90 
(90 % confidence range) times the cosmic value. 
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Fig. 4. — The spectra of GX 301—2 from 5 keV to 10 keV for the first and second halves of the 
PP-phase observation. The solid lines show the Monte Carlo models inferred from the data. 
In addition to the iron Ka emission line profile, the continuum shape is also successfully 
reproduced with the parameters inferred from the Compton shoulder profile. The values for 
A^h are 12.0 x 10 23 cm~ 2 and 8.5 x 10 23 cm~ 2 for the first and second halves, respectively. 
A photon power-law index of 1.0 is assumed. 



